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ABSTRACT 

The Laser Interferometer Space Antenna (LISA) will provide the largest observational sample of (interacting) 
double white dwarf binaries, whose evolution is driven by radiation reaction and other effects, such as tides 
and mass transfer. We show that, depending on the actual physical parameters of a source, LISA will be able 
to provide very different quality of information: for some systems LISA can test unambiguously the physical 
processes driving the binary evolution, for others it can simply detect a binary without allowing us to untangle 
the source parameters and therefore shed light on the physics at work. We also highlight that simultaneous 
surveys with GAIA and/or optical telescopes that are and will become available can radically improve the 
quality of the information that can be obtained. 

Subject headings: gravitational waves - instrumentation: detectors - instrumentation: interferometers - meth- 
ods: data analysis - binaries (including multiple): close - white dwarfs 



1. INTRODUCTION 

White dwarf (WD) binary systems are the most abundant 
class of compact binaries in the Galaxy. Double white dwarf 
binaries (DWDs) come in two distinct flavours: detached bi- 
naries (WD-WD), that are driven to shorter and shorter orbital 
periods primarily by loss of angular momentum through grav- 
itational wave radiation and interacting white dwarf binaries 
(AM CVn systems), where matter is being transferred from 
the lower-mass white dwarf (the donor star) to the more mas- 
sive accretor. These two classes can be unified in one picture, 
in which detached systems evolve to shorter periods with tidal 
effects becoming progressively more important, until mass 
transfer ensues. Depending on the mass ratio of the WD-WD, 
the system will either merge (when the mass ratio is larger 
than about 2/3) or reach a "turning-point" after which the pe- 
riod evolution is reversed and the binary evolves to longer pe- 
riods, with ever decreasing mass ratio (e.g. Natheret al. 1981, 
Tutukov & Yungelson 1996, Nelemans et al. 2001a, Marsh et 
al., 2004). 

These binaries are important laboratories for some funda- 
mental questions in astronomy: binary evolution, the structure 
of WDs, tidal effects as well as the nature of the progenitors 
of type la supernovae (e.g. Tout 2005). Despite their large 
population, the present sample of WD binary systems is fairly 
limited; WDs are electromagnetically faint and spectroscopic 
surveys suffer from strong observational biases, in particular 
against systems characterised by short orbital periods, which 
disappear quickly from the observable sample through merg- 
ing. The ESO SPY survey (Napiwotzki et al. 2004) has un- 
dertaken an extensive spectroscopic survey to search for WD- 
WD systems with the primary aim of testing the Supernova la 
scenario. This has led to the detection of about ~ 100 WD- 
WD systems, which represents a « 10-fold increase of the 
number of known systems with respect to the pre-SPY era 
(Marsh 2000). The 17 known AM CVn systems are all found 
as by-products, mainly from quasar surveys (see Warner 1995 
and Nelemans 2005 for reviews). 

The Laser Interferometer Space Antenna (LISA) (Bender 



et al, 1998), an ES A/NASA gravitational wave observatory 
with launch date 2013+ is expected to dramatically increase 
the number statistics of observed white dwarf binaries by a 
factor Zt 10 by observing the gravitational radiation emitted 
by these systems. LISA will provide the largest observational 
sample of these objects, with 10 3 - 10 4 systems detected dur- 
ing the mission lifetime (e.g. Nelemans et al, 2001b); in addi- 
tion, such survey is not hampered by the same biases of those 
carried out in the electromagnetic band: LISA will probe the 
whole galaxy (and possibly beyond to reach the Small and 
Large Magellanic Clouds) and will be very effective in iden- 
tifying short period ( ^ 30 min) binaries. A handful of in- 
teracting binaries known from previous surveys are so close 
and with sufficiently short periods that gravitational waves 
(GWs) emitted by them can be confidently detected within 
a few weeks-to-months of the operation of the instrument. In- 
deed DWDs are guaranteed sources for LISA. It is also ex- 
pected that a fair fraction of the LISA detected systems will 
be observable with X-ray, optical and possibly astrometric in- 
struments (such as GAIA), largely increasing the amount of 
information that can be extracted (e.g. Nelemans et al. 2004, 
Cooray et al, 2004). 

Despite the crucial role of DWDs in LISA observa- 
tions, studies of GW astronomy have been devoted to either 
monochromatic detached WD-WD or WD binaries where the 
drift in frequency is assumed linear and "clean" (e.g. Taka- 
hashi and Seto, 2002; Seto 2002), in the sense that it is ex- 
actly described by the evolution of two point-masses in vac- 
uum. In fact, the astrophysical scenario is radically differ- 
ent: by the time the frequency evolution becomes important, 
tidal and other effects play such an important role that on the 
one hand the frequency evolution offers a unique opportunity 
to studying these effects, and on the other hand astronomical 
observations that have been suggested possible with LISA - 
such as the direct measure of the source distance D and chirp 
mass Ai - will be much more difficult, because of the chal- 
lenge of disentangling in an unambiguous way several physi- 
cal effects. 

This raises two key questions that are addressed for the 
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first time in this Letter: (i) Can LISA test unambiguously the 
physical process(es) driving the evolution of a WD binary? 
As a consequence, (ii) What is the information that can be 
gained from the data, if one drops the a priori (and unrealis- 
tic) assumption that radiation reaction is the only mechanism 
at work? The answers to these questions depend on a wide 
range of factors that we discuss, in particular the source period 
and mass and whether astrometric and/or optical observations 
are available to complement LISA. 

2. ASTRONOMY WITH LISA 

WD binaries enter the LISA observational window (0.1 
mHz - 100 mHz) when they reach a period P w 5 hr; they 
secularly evolve through radiation reaction across the LISA 
band and are so numerous to create a stochastic foreground 
that dominates the LISA observational window up to m 3 mHz 
(e.g. Hils et al, 1990). When a binary has reached a gravita- 
tional wave frequency / = 2/P rts a few mHz (and can there- 
fore be resolved individually), the recorded signal, over a typ- 
ical observation time T = 1 yr, shows an intrinsic frequency 
evolution. Depending (primarily) on the mass ratio, the fate 
of the binary can be dramatically different. Comparable mass 
stars continue their inspiral until they become unstable and 
merge in the frequency band 10-100 mHz (depending on their 
mass). If the lower mass companion is sufficiently light, sta- 
ble mass transfer can start. The orbit eventually stalls - the 
frequency evolution as a function of time reaches a "turning 
point" - and then the period increases with secular evolution 
dominated by mass transfer: the binary becomes observable 
as an AM CVn system. 

LISA is an all-sky observatory; it monitors DWDs (and 
other sources) providing two independent data sets synthe- 
sised using the technique known as time-delay interferome- 
try (TDI), see Dhurandhar and Tinto (2005) for a recent re- 
view. Each data set (that we label with a = 1,2) can be for- 
mally represented as s (a \t) = h (a \t;\) + n (a) {t), where h (a \t ; A) 
is the GW strain, A the vector of the unknown signal param- 
eters and n (a) (f) the noise. Here we model the gravitational 
waveform at the lowest Newtonian quadrupole order - in fact 
post-Newtonian corrections to the phase contribute <C 1 wave 
cycle for multi-year observations and harmonics at other mul- 
tiples of 1 /P are a factor <, 10 3 smaller than the leading order 
Newtonian quadrupole (Blanchet et al, 1996) - and the noise 
according to the analytical fit by Barack and Cutler (2004), 
ignoring the possible contribution from unresolved extreme 
mass ratio in-spirals. We model the detector output using the 
so-called rigid adiabatic approximation (Rubbo et al, 2004), 
so that the signal registered at LISA outputs can be repre- 
sented as (Vecchio & Wickham, 2004) 

4 

h (a \t)=A Y / A ( :\t)cos X ^(t), (1) 
n=l 

where 

A// 5 / 3 

Ao = 2(Vo) 2/3 -^-; (2) 

Ajf^jN, t,,f(t)] and cosxl a) I>;N,L,/(f)] contain informa- 
tion about the intrinsic frequency evolution /(f) - here fo is 
the frequency at the beginning of the observations - and the 
source position in the sky N and the orientation of its orbital 
plane L; explicit expressions for h {a \t) can be found e.g. in 
Vecchio & Wickham (2004). 



The intrinsic frequency of radiation from a DWD slowly 
changes during the typical LISA observation time. From a 
phenomenological point of view, we can therefore Taylor ex- 
pand it around an arbitrary fixed time (say the beginning of 
the observations fo = 0), so that the phase of the signal can be 
written as: 

JL t (k+V) d k f 
^ = ^2^ + 2^ — ^, (3) 

where we define with f {k) = d k f/dt k (k = 1, ...,K) the "spin- 
down parameters" - notice that despite the terminology f {k) 
could be either positive or negative. The number of spin-down 
parameters that is retained into the phase model depends on 
the actual source parameters and is determined by the con- 
dition that the integrated value of the phase over the whole 
observation time differs from the actual one by -C 1 . If both 
fo and fo are observable (K = 2), then one can measure the 
so-called braking index n = fofo/fo, which provides a pow- 
erful tool to discriminate the physical mechanism driving the 
orbital evolution; if only radiation reaction is at work, then 
n=ll/3. 

The signal h {a \t ; A) depends therefore on a vector A of 7+K 
unknown parameters: the amplitude Aq, four angles that iden- 
tify N and L, the frequency at the beginning of the observation 
fo, an overall arbitrary phase <f>o and K spin-down parame- 
ters. We can then estimate the expected mean-square-error 
((AA-0 2 ) (j labels the components of A) which is associated 
to the LISA measurement of each parameter X j by computing 
the variance-covariance matrix (Vecchio & Wickham, 2004); 
the diagonal elements of the above matrix provide (in the limit 
of high signal-to-noise ratio p) a tight lower limit to ((AA J ) 2 } 
(see, e.g. Nicholson and Vecchio, 1998). 

A thorough exploration of the parameter space characteris- 
ing DWDs goes beyond the scope of this paper. Here we con- 
centrate on two fiducial sources that are representative of the 
astrophysical scenarios considered here (Marsh et al, 2004): 
(i) An "unstable" binary with m\ = ni2 = O.6M , that never 
reaches a turning point and merges at 37.7 mHz; in order to 
faithfully describe the phase evolution one needs to include fo 
for fo > 0.3 mHz (1 mHz) and both f and f for / > 8 .8 mHz 
(21.7 mHz) for an integration time T = 5 yr (1 yr); (ii) A "sta- 
ble" binary with individual progenitor masses m\ = 0.4M Q 
and ni2 = 0.2M Q , whose frequency increases until it reaches 
a turning point at / = 9.7 mHz and then enters the AM CVn 
phase where the frequency decreases; the phase evolution is 
described by only the first spin-down parameter fo in the fre- 
quency range 0.5 mHz (1.2 mHz) < fo < 9.7 mHz (in the in- 
spiral phase) and 0.8 (2.6 mHz) < fo < 9.7 mHz (in the AM 
CVn phase) for T = 5 yr (1 yr); elsewhere the source is de- 
tected as monochromatic. In our calculations we take tidal ef- 
fects and mass transfer explicitly into account. However, we 
do not consider short timescale variations that are observed in 
many binaries and which might have a (strong) influence on 
the instantaneous value of spin-down parameters (e.g. Marsh 
& Nelemans, 2005). 

In order to quantify how accurately LISA can measure the 
parameters of interacting WDs, we have carried out Monte 
Carlo simulations over a population of 10 4 sources charac- 
terised by the same physical parameters -Aq, fo, fo and fo, if 
relevant - but with random values of N and L. All the sources 
are placed at the same distance D = 10 kpc (which in turn 
sets Aq). We simulate a population of unstable binaries at 
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TABLE 1 

SlGNAL-TO-NOISE RATIO p AND THE STATISTICAL ERRORS 
ASSOCIATED TO THE MEASUREMENT OF THE AMPLITUDE A, THE 
POSITION OF THE SOURCE IN THE SKY AQ N , THE INCLINATION COS i, 
THE SPIN-DOWN PARAMETERS, f AND fij, AND THE BRAKING INDEX n 
IN OBSERVATION OF AN UNSTABLE BINARY AT 30 MHZ (TOP) AND 20 
MHZ (MIDDLE) AND A STABLE AM CVN SYSTEM AT 8 MHZ (BOTTOM). 

WE LIST THE MEDIAN AND THE VALUE CORRESPONDING TO A 
CUMULATIVE PROBABILITY OF 90% OF THE SEPARATE DISTRIBUTIONS 
FOR EACH PARAMETER. IN EACH CASE WE SIMULATE 10 4 SOURCES AT 
DISTANCE D = 10 KPC WITH RANDOM SKY POSITION AND ORIENTATION 
OF THE ORBITAL PLANE. THE OBSERVATION TIME IS T = 1 YR AND 5 YR. 





T 


P 


AA/A 


AQ N 


A cost 


A/0//0 


A/0//0 


An 




(yr) 






(deg 2 ) 










median 


1 


23 


0.13 


0.09 


0.31 


1.3 x 10~ 3 






90% 


1 


14 


0.30 


0.29 


0.68 


2.3 x l(r 3 






median 


5 


52 


0.06 


0.004 


0.2 


4.1 x 10~ 5 


8.1 x ur 3 


0.03 


90% 


5 


32 


0.14 


0.01 


0.4 


6.5 x i(r 5 


1.3 x ur 2 


0.05 


median 


1 


33 


0.09 


0.12 


0.26 


4.4 x 10~ 3 






90% 


1 


20 


0.20 


0.31 


0.56 


7.4 x 10~ 3 






median 


5 


74 


0.04 


0.005 


0.17 


1.3 x 10~ 3 


7.6 x ur 2 


0.28 


90% 


5 


45 


0.09 


0.01 


0.37 


2.1 x 10~ 3 


1.2 x ur 1 


0.45 


median 


1 


8.1 


0.35 


11.75 


0.52 


1.6 






90%, 


1 


5.0 


0.82 


30.93 


1.10 


2.8 






median 


5 


18 


0.15 


0.43 


0.34 


8.0 x 10~ 3 






90% 


5 


11 


0.40 


1.08 


0.74 


1.3 x 10~ 2 







fo = 20 mHz and fo = 30 mHz, respectively, where the orbital 
separation is « 5 star radii. These frequencies - chosen to 
provide a "best case scenario" from the point of view of in- 
formation extraction - are at the very high end of the range 
expected for Galactic DWDs (Nelemans et al, 2001b). We 
also simulate a population of stable binaries at /o = 8 mHz in 
the AM CVn phase, where mi = 0.43 M , m 2 = O.17M and 
fo = -1.2 x 10~ 15 s~ 2 . We compute the variance-covariance 
matrix and summarise the results of the distributions of the 
expected mean-square-errors associated to the LISA measure- 
ments in Table 1 . 

We consider now what are the consequences for astronomy 
and concentrate on the specific results for T = 5 yr; the reader 
can easily re-scale the conclusion for T = 1 yr, using the re- 
sults provided in Table 1. About a third of the DWDs de- 
tectable by LISA are also expected to be observed as eclips- 
ing binaries (Cooray et al, 2004) yielding a direct measure of 
R\/(a cosi), Ri/ia cost) and R1/R2, where R12 are the radii 
of the two stars, a is the orbital separation and cos 1 = N • L 
the orbital inclination. In addition, about half of the DWDs 
might be detected by GAIA; according to Cooray et al. (2004) 
the distribution of the V-magnitude of galactic DWDs peak 
around 20 which corresponds to the magnitude limit up to 
which GAIA can do high precision observations 1 . The astro- 
metric accuracy of the instrument is 200 /ias for V = 20, 10/ias 
for V = 15, and 4/ias for V < 12 which implies that GAIA can 
detect DWDs up to a few kpc and measure their distance with 
an error in the range « 1%- 10%. The possibility of observ- 
ing DWDs with LISA, GAIA and/or spectroscopically may 
be crucial in disclosing the physics at work in these systems 
by breaking the degeneracy of a large number of parameters, 
as we now discuss in more detail. 

We focus first on the unstable DWD for which both fo and 

1 However, the estimate by Cooray et al. (2004) is based on very optimistic 
calculations of tidal heating by Iben et al. (1998) 



fo are observable. This is crucial, because one can measure 
the braking index with an absolute error An 0.1, and LISA 
will be able to test directly the physical processes that drive 
the binary evolution. If observations are consistent with the 
assumption that the evolution is driven by radiation reaction 
only, then from fo = (96/5) 7r 8 ' 3 .M 5 / 3 / 11 / 3 one can measure 
directly the chirp mass of the binary with a statistical error 
{(AM/M) 2 } 1/2 ~ (3/5)((A/ //o) 2 ) 1/2 ~ 10" 4 and, through 
Eq. the distance with ((AD/D) 2 } [ / 2 ~ ((AA/A) 2 ) 1 / 2 « 0.1 
(in both cases we have ignored the errors on fo and the corre- 
lations that are negligible). If a LISA source is also observed 
as an eclipsing binary, then one is able to measure R\ and R2 
(cosi is estimated with a good accuracy w 10% by LISA) and 
place an upper-limit on the total mass of the system: these 
combined observations represent a unique opportunity to dis- 
close the physics of compact objects. 

However, our theoretical understanding of DWDs, despite 
limited, suggests that by the time a source emits in the mHz 
region the orbital evolution is not clean, which, for the par- 
ticular source we are considering now, will be unveiled by 
LISA through a direct measure of n not consistent with 11/3. 
In this case fo and fo carry the signature not only of A4 but 
also of the parameters that control tidal effects. If only LISA 
data are available, then neither A4 nor D can be inferred (only 
loose limits can be placed): this is in strong contrast with 
what has been assumed so far, based on the naive hypothesis 
that hydrodynamical effects can be neglected. Nonetheless, 
the situation is radically different if the LISA source is also 
observed by GAIA which provides an independent measure 
of D: from it and the amplitude of the gravitational signal, 

Eq. one can estimate M. with an error ( ( AM / M) ) 1//2 = 

(3/5)\/ '{(AD/D) 2 )+ ((AA/A) 2 ) w 0.1, if one assumes that 

((AD/D) 2 ) 1 / 2 « ((AA/A) 2 ) 1 / 2 w 0.1; for some systems one 
could achieve an accuracy as good as w 0.01. In other words, 
GAIA re-instates the full power of LISA observations by 
breaking the degeneracy between M. and D. From M, fo and 
fo one can learn new physics that we have not been able to 
explore directly so far. Of course the most exciting situation 
is if the LISA and GAIA source is also identified optically as 
an eclipsing binary, which provides additional information on 
the stars size and total mass. 

We turn now attention on the observations of the stable AM 
CVn progenitor. The key difference from the scenario de- 
scribed above is that LISA is now able to detect only fo; more- 
over the signal-to-noise is much smaller, which degrades the 
quality of astronomy in general, and the size of the LISA er- 
ror box in the sky in particular: the association of a source 
observed by LISA with stars detected by other surveys be- 
comes much more challenging. If LISA measures a value 
fo < 0, then this provides a strong evidence that the DWD 
evolution is driven by mass transfer; in this case however, D 
can not be inferred because fo depends on A4 and other phys- 
ical parameters; the chirp mass is equally unmeasurable. On 
the other hand, if fo > 0, then one is not able to make any 
statement about whether the system is "clean" or not. An in- 
dependent determination of the distance (e.g. through GAIA) 
becomes again essential: if this is possible, one can compare 
the value of fo predicted by evolution driven purely by radia- 
tion reaction with the value obtained with LISA. For AM CVn 
the ability of testing in detail models of orbital evolution and 
mass transfer is somewhat reduced because only one parame- 
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ter (/o) is directly measured; from this point of view, a mission 
lasting up to 10 years (which is not a mission requirement at 
present but is considered feasible with the instrumentation on 
board) is crucial, as fo becomes observable for an increasingly 
larger number of sources. It is also worth emphasising that if 
an AM CVn is observed optically, then one might measure 
the super-hump period. Based on extrapolation of an empir- 
ical relation, this yields the mass ratio of the binary with an 
accuracy of a few tens of percent (e.g. Espaillat et al. 2005). 
From A4 (provided by LISA and GAIA) and m\/m.2, one can 
determine the individual masses of the stars and, combined 
with /o, the mass transfer rate. 

3. CONCLUSION 

Interacting WD binaries are unique laboratories for test- 
ing theories about binary evolution, the structure of white 
dwarfs, tidal effects, mass transfer and supernova-la progen- 
itors. LISA opens a radically new opportunity to carry out 
these studies by providing a much larger (and unbiased) sam- 
ple of these stars and information on physical parameters that 
are hardly accessible by other surveys. However the challenge 
that LISA astronomers will face is to untangle a large number 
of physical parameters from a few phenomenological observ- 
ables. 

We have shown that if the first two time derivatives of the 
frequency, fo and fo, are measured, then LISA alone will be 
able to test accurately whether or not a system evolves cleanly 
due to radiation reaction. If the evolution is driven by the 
competition of several processes (which we expect in a large 
number of cases), then the phase evolution monitored by the 



instrument allows us to constrain models of evolution of WD 
binaries. If only fo is observed, then LISA provides no infor- 
mation on the problem at hand. If the evolution is not "clean" 
and only LISA data are available, one is not able to infer the 
chirp mass and distance of the binary: this is in strong contrast 
with what has been assumed so far based on the simplistic as- 
sumption that tidal effects and mass transfer can be ignored. 
However, if a LISA source is also observed by GAIA the sit- 
uation changes dramatically: the combined observations pro- 
vide both M. and D with an accuracy » 10%. Moreover com- 
bining LISA and GAIA observations, it would still be possi- 
ble to test whether the binary evolution is clean even if LISA 
can measure only fo- In addition, a fair fraction of DWDs de- 
tected with LISA are observable as eclipsing binaries, from 
which one can determine the stars individual radii: this opens 
the possibility of extracting an unprecedented wealth of infor- 
mation about the physics of these compact objects in the most 
extreme physical conditions. 

Our study calls for a more detailed investigation of this ob- 
servational scenario over the full parameter space (such work 
is presently in progress and will be reported elsewhere) in- 
cluding the exploitation of the full power of the TDI observ- 
ables (e.g. Dhurandhar et al, 2002, Prince et al, 2002) and 
detailed theoretical modelling of interacting white dwarfs in 
realistic conditions in order to be able to link reliably the phe- 
nomenological observational parameters to those that under- 
pin the key physics. 
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